'AClUry FORM 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


, / • . 


Technical Report 32Tt$9 




Determination of Aerodynamic Damping Coefficients 
From Wind-Tunnel Free- Flight Trajectories 
of Non-Axisyrnmetric Bodies 


R.H. Pris/in 
M.B. Wilson 



JET PROPULSION LABORATORY 

CALIFORNIA INSTITUTE OF TECHNOLOGY 
PASADENA, CALIFORNIA 


October 15, 1967 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


Technical Report 32-1159 

Determination of Aerodynamic Damping Coefficients 
From Wind-Tunnel Free-Fiight Trajectories 
of Non-Axisymmetric Bodies 

R.H. Pris/in 
M.B. Wilson 


Approved by* 


E. A. Laumann, Manager 
Aerodynamic Facilities Section 


JET PROPULSION LABORATORY 

CALIFORNIA INSTITUTE OF TECHNOLOGY 
PASADENA, CALIFORNIA 


October 15, 1967 


TECHNICAL REPORT 32-1159 


Copyright © 1 967 
Jet Propulsion Laboratory 
California Institute of Technology 

Prepared Under Contract No. NAS 7-100 
National Aeronautics & Space Administration 



Contents 


I. Infroduction 1 

A. Preliminary Remarks 1 

B. Some Elementary Equations of Statics 2 

C. Equations of Planar Motion 2 

D. Scope of the Report 7 

II. The Energy Integral Relation and Its Application 8 

A. An Integral of Planar Motion 8 

B. Offset Linear Pitching Moment: C m ($) = C„ J() + C,„ a 0 10 

C. Bi-Linear Pitching Moment 11 

D. Bi-Cubic Pitching Moment 12 

E. Full Cubic Pitching Moment 16 

III. Equivalence Relations Between Moment Curves 21 

A. Equivalence Criteria 21 

B. Offset Linear and Bi-Linear Pitching Moments 22 

C. Bi-Cubic Pitching Moment 22 

D. Cubic Pitching Moment: C w (0) = C m j9 + 3/2 C W] O' 1 + 2 C m> 0* 25 

IV. Application and Accuracies 26 

A. Approach 26 

B. Sample Case 1 26 

C. Sample Case 2 29 

Appendix: Iterative Formulas for the Solution of f (X) ~ 0 33 

Nomenclature 33 

References 34 

Tables 

1 . Summary of distance periods over the V 4 cycfes 14 

2. G 3 vs a 2 and k 2 . 19 

3. Comparison of results using various pitching moment approximations 

for sample case 1 29 

4. Comparison of results using various pitching moment approximations 

for sample case 2 31 

JPL TECHNICAL REPORT 32- FI 59 iii 


Contents (contd) 


Figures 

1. Definition sketches for force and moment conventions 3 

2. Schematic representations of body shapes and corresponding 

static aerodynamics 4 

3. Data reduction coordinate system 5 

4. Comparison of motion in coordinate and phase planes 8 

5. Complete elliptic integral of the first kind vs elliptic modulus 15 

6 . Errors introduced by using a linear moment to approximate a 

cubic moment 15 

7. The functions Gi and G 2 vs elliptic modulus 17 

8. Static aerodynamics for sample case 1 26 

9. Roots of the potential energy equation for sample case 1 27 

10. Elliptic integral parameters for sample case 1 27 

11. Trajectory characteristics for sample case 1 27 

12. Integral of 6' for sample case 1 27 

13. Lift and drag terms for sample case 1 28 

14. Six-degree-of-freedom angular history, sample case 1 28 

15. Equivalent bi-cubic pitching moment, sample case 1 29 

16. Static aerodynamics for sample case 2 29 

17. Trajectory characteristics for sample case 2 30 

18. Negative vs positive peaks for sample case 2 30 

19. Amplitude decay history for sample case 2 30 

20. Lift and drag terms for sample case 2 30 

21. Equivalent bi-cubic pitching moment, sample case 2 31 

22. Equivalent cubic pitching moment, sample case 2 31 


iv 


JPL TECHNICAL REPORT 32-1159 



Abstract 


General formulas for estimating the aerodynamic damping coefficient (C m<? + C m 
from planar free-flight oscillation histories are presented. Using Jaffe's energy- 
integral equation (a phase plane integral of the motion), systematic approximations 
are introduced that lead to the derivation of accurate results for highly nonlinear 
and asymmetric static aerodynamics. Particular attention has been given to making 
the results valid, even for large initial amplitudes and for extreme asymmetries. 
Two examples are presented to demonstrate the accuracies involved with the 
various approximations and to illustrate the computational procedure needed for 
actual data reduction. 
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Determination of Aerodynamic Damping Coefficients From 
Wind-Tunnel Free-Flight Trajectories of 
Non-Axisymmetric Bodies 


I. Introduction 

A. Preliminary Remarks 

Free-flight testing techniques have become increasingly 
useful in the wind-tunnel study of high-speed aerodynamic 
characteristics. During the past 4 yr, the Jet Propulsion 
Laboratory (JPL) has put considerable effort into the 
development of these techniques for use in conventional 
wind tunnels. Reference 1 provides a good survey of the 
work accomplished to date. 

The primary value of free-flight testing lies in the 
absence of model support interference. Furthermore, the 
technique allows an experiment to be conducted in a 
dynamic environment similar to actual vehicle flight, 
rather than in the standard static mode. Therefore, one 
very important application of such testing techniques 
is the determination of aerodynamic force and moment 
coefficients and, in particular, obtaining the aerodynamic 
damping moment coefficient (C Mq + C m& ) 9 more often 


called the dynamic stability derivative. Studies at JPL 
and elsewhere have generated successful methods of 
deducing this damping coefficient from the planar oscilla- 
tory history provided by free-flight test data, (e.g., Refs. 2, 
3, and 4). Thus far, the studies have dealt exclusively with 
axisymmetric configurations that have their centers-of- 
gravity on the centerline axis of the bodies. The cor- 
responding symmetry in the vehicle aerodynamics leads to 
certain simplifications of trajectory analysis and the result- 
ing data reduction. However, this is a restriction that pre- 
cludes the possibility of analyzing data corresponding to 
asymmetric configurations, such as lifting bodies, as well 
as axisymmetric shapes with offset centers of gravity. 

Tt! onfi/^iT\otiArk rvf fnfnra *■ p ft*pp flirflif fprfirwr 

iii v/i iuiuiv imviV/Oi m liiC iiC/C~iiigiii iC5uii£ 

of such configurations, this report discusses some initial 
studies aimed at the problem of obtaining dynamic sta- 
bility derivatives from the planar free-flight trajectories 
of non-axisymmetric vehicles. We make the restriction of 
planar motion and, in so doing, confine the report to bodies 
with a plane of symmetry, but not necessarily an axis of 
symmetry. 
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B. Some Elementary Equations of Statics 

To determine dynamic coefficients from free-flight tra- 
jectory data it is necessary to have complete knowledge 
of the vehicle's static aerodynamics. In general, the static 
coefficients can be obtained from the free-flight trajectory 
prior to dynamic analysis (Ref. 2). In this report, it will be 
assumed that these coefficients are known. Because the 
static aerodynamics can be given in several forms, it is 
useful to specify which notations have been used. Consider 
an arbitrarily shaped body (with a plane of symmetry) at 
static equilibrium in a uniform flow, V M . Because the 
motion is confined to one plane, say the vertical plane, 
only two directions are needed to specify the forces acting 
on the body. Define the axial and normal forces, F A and 
F v , as shown in Fig. la. Lift L and drag D are defined in 
the usual fashion. Resolving these forces, we get the fol- 
lowing two sets of equations relating the aerodynamic 
coefficients C h , C [h C v , and C t : 

C L = C, cos a — C A sin « 

C„ = Cv sin a + C, cos a 

C v = C h cos a -b Cd sin « 

C t = — C /y sin<* + C D cos a ) 

The same relations apply for an axisymmetric shape with 
an offset center of gravity (eg) (Fig. lb). In either case, the 
reference axis is taken through the eg so that the dynamic 
equations remain as simple as possible. 

Summation of the moments about the vehicle's eg gives 



angle-of-attack, «. As it is not desired to put restrictions on 
the vehicle oscillation amplitude, it is anticipated that the 
cp location will also vary with <*, i.e., ( l cp — l cg ), as well 
as e z , may be functions of a. Nevertheless, the Eqs. (1-4) 
may always be applied at any fixed a to relate the static 
coefficients. 


Figure 2 shows the general types of static aerodynamics 
considered in this report. Note that for an axisymmetric 
body with an offset eg, the offset has no effect on the lift 
and drag curves. However, the moment curve is altered so 
that the axisymmetric methods of Refs. 2, 3, and 4 do not 
apply. No restrictions on the linearity of the coefficients 
are imposed. 


C. Equations of Planar Motion 

The coordinate system used for the analysis is shown 
in Fig. 3, where x m — z m is an inertial axis system, fixed 
with respect to the wind tunnel. It is with reference to this 
axis system that all measurements of the model's positions 
are made. However, for analysis purposes, distances must 
be equivalent to those traveled by a corresponding vehicle 
moving through still air. The measured value in the z m 
direction is correct, but the x m distance is incorrect because 
of the wind-tunnel’s free-stream velocity. Therefore, the 
X-Z axis system is defined as fixed with respect to the mov- 
ing media. Prior to data reduction, the following trans- 
formation equations must be applied (Fig. 3): 


X — x m + Vo Q t 
X = x m + Voo 


M v9 = F, (e 7 ) - Fv (In - u (3) 

or in coefficient form 

C. = C,(f)-C,(^) (4) 

where e z is the center of pressure (cp) offset normal to the 
reference axis, d is the reference length, and l ( . t , — l r</ is the 
distance between the cp and eg along the reference axis. 

For any particular vehicle configuration, the static coeffi- 
cients C L , C/j, and C m are, themselves, functions of the 


X = x„, (5) 

Because of the aerodynamic lift force, a model velocity 
in the Z-direction is induced; z„, ~ Z. The actual vehicle 
velocity V in the X-Z coordinate system is the vector sum 
of the three velocities Voo, x,„, z,„ and, in general, will not 
be parallel to the X-axis. The model angle-of-attack, «, 
defined as the angle between the velocity vector and the 
model reference axis, cannot be measured directly from 
wind-tunnel data. Therefore, the equation of angular 
motion will be written in terms of the more convenient 
angle 0 (Fig. 3). Portions of the ensuing derivation of 
this equation of angular motion are substantially equiva- 
lent to the derivation of Ref. 4. However, the entire 
development is included here for purposes of clarity and 
completeness. 
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LIFT CURVE 1 LIFT CURVE 

(EXACTLY ANTISYMMETRIC) 



(EXACTLY SYMMETRIC) 



(ANTISYMMETRIC CURVE 
DISPLACED BY C mQ 


Fig. 2. Schematic representations of body shapes and corresponding static aerodynamics 
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Fig. 3. Data reduction coordinate system 
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Substituting in Eq. (9), 


9" = + (C mq + C m .) e’d - (O' - a') C m .d + ^C D e'^ ( 12 ) 

The static pitching moment C m is a function of From expanding the function C,„ (a) about 9 in a Taylor series and 
ignoring the higher order powers of (9 — a). 


C„,(a) — C,„(ff ) + (a ~ 9) ^0 ^ 


(13) 


From differentiating Eq. (11) and using Eq. (7), 


Z = X(9-a)+X(6 


a) — — 


pAX- 

2m 


C n (9 - a) + X-{9' - a') 


(14) 


Combining Eqs. (8) and (14), 


a' -9'=- + 4- + ■£— C i) (a - 9) 

2m x- 2m ’ 


(15) 


The third term on the right side is about two orders of magnitude smaller than the first and it will be ignored. Integrat- 
ing Eq. (15) with respect to X, 


(a ~ 0) — (a ” 0)o + 



pA ^ 

“ 2 ^ C/ ' + 



(16) 


where — 0 O is the angle corresponding to X = 0. 


The lift coefficient C L , which appears in Eq. (16), and the drag coefficient C D , which appears in Eq. (12), are also 
functions of a. However, they comprise second-order effects in the moment equation and, because 0 is about equal to a, 
they can be assumed to be functions of 6 directly with a negligible error resulting. Also, normally | O' | > > | 0' — a\ 
and in general | C niq \ >> |C m -|; therefore, the third term on the right side of Eq. (12) is negligible. By control of the 
initial conditions, i.e., Z 0 = 0, we can have (a — 0) o — 0. Finally, for a wind-tunnel free-flight application, gravity terms 
have no appreciable effect. Incorporating these facts into a combination of Eqs. (12) and (16) yields the following 
integral-differential equation of angular motion: 


(9" 


pAd pAd / pA\ dCJfi) [" CjW 

2i C m {0)+ 2I ^ 2m J d() J' e , 




pAd- 

2/" 


(C„, + C m .) 9 ’ = 0 


(17) 


‘The static aerodynamic coefficients will, in general, also vary with the Mach number. However, because of the extreme dominance of 
the V„ term in the model velocity (supersonically about 98%), the Mach number is assumed to remain constant for the flight. 
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The problem, then, consists of solving the nonlinear 
integral-differential, Eq. (17), and inverting the solution 
to obtain (C Wq -f C m .) as a function of 0 and X. The accu- 
racy of the final expression for the damping derivative 
depends on the validity of the assumptions made both in 
the development and solution of Eq. (17). It is, therefore, 
instructive to summarize the assumptions used to this 
point and, to mention some common restrictions usually 
placed on this type of problem but which have not been 
or will not be imposed in this report. 

With the exception of dropping the gravity term, there 
have been only two assumptions made. The first is that 
(0 — «) is a very small angle. The effects of this small 
angle have been ignored in second-order (and higher) 
terms, and handled in an approximate manner in the pri- 
mary term. The second assumption is that the static pitch- 
ing moment is the main moment acting, and all other 
terms of the moment equation have second-order effects. 
For a wind-tunnel application, these two suppositions do 
not result in any significant errors, even when integrations 
are performed over long time periods. Therefore, Eq. (17) 
is nearly an exact expression for planar angular motion. 


suiting from the energy-integral approach. Therefore, the 
work presented in this report consists primarily of appli- 
cations of the energy integral equations to various sets of 
nonlinear aerodynamics. 

To solve Eq. (17), it is necessary to have explicit func- 
tional relationships between the static aerodynamic coeffi- 
cients and the angle 0. Because any continuous function 
can be represented as a power series expansion, it is natu- 
ral to consider many-termed polynomials for C D> C L , and 
C m . Lift and drag are of second-order influence on the 
motion and, therefore, such representations offer no par- 
ticular problems, as approximations may be made in the 
evaluation of resulting integrals without significant loss 
of accuracy. However, the same freedom of approxima- 
tion is not available when considering the primary term: 
the pitching moment. Consequently, the main limitation 
on nonlinearities and the corresponding generality of the 
solution must be imposed on the pitching moment. 

Throughout the remainder of this report, lift and drag 
coefficients will be considered with the following func- 
tional representations: 


No mention has been made of the functional depen- 
dence of (C niq + C m >). The dynamic stability derivative is 
usually regarded as an effective constant coefficient over a 
fixed time period (e.g., the flight time, one-half the flight 
time, one oscillation cycle, etc.). This is not really a lim- 
itation because, by testing at various amplitudes (or fre- 
quencies, etc.), an effective coefficient vs amplitude curve 
can be obtained and, if desired, from this the damping 
derivative (as a function of instantaneous angular dis- 
placement) may be determined (Ref. 5). No restrictions of 
linearity of static aerodynamics or small angular excur- 
sions have been made in the derivation of Eq. (17). 


c, (0) = c Lo + c La o + c L p + C h Jt' | 

C/,(0) = C Do + C Dl 0 + Cd 2 0 2 + Cd 3 0 3 + C D J 4 j 


(18) 


Additional terms for these series expansions could be han- 
dled without too great of an increase in complexity. How- 
ever, these particular curves are sufficiently general to be 
representative of a large class of shapes, and more terms 
would not be necessary. The following four alternate 
pitching moment curves are considered in detail: 

(1) Case 1, offset linear moment: 


D. Scope of The Report 

The primary motivation for the present study was to 
develop analytical tools capable of handling cases of non- 
axisymmetric bodies or axisymmetric bodies with a eg 
offset. The established energy-integral approach used in 
Refs. 2, 3, and 4 was regarded as one possible method, 
but not necessarily the only or best approach. Consider- 
able effort was expended into investigating alternate pro- 
cedures and fresh viewpoints for the solution of the 
differential Eq. (17). Examples are small perturbation 
techniques (Ref. 6) and transform methods that simplify 
the original differential equation (Ref. 7). However, all 
other approaches investigated led to more serious approx- 
imations and more complicated solutions than those re- 


c m (0) = Cm 0 + C m<x 0 

(2) Case 2, bi-linear moment: 


c m (o)=-c ma _e 9 o^o ^ 

c m (0) = c ma j, e^o j 

(3) Case 3, bi-cubic moment: 

C m (6) = Cm a 6 + 2C m2 0 3 , 6^0 

Cm (9) = Cm a 6 + 2Cm 1+ 6 3 , 6^0 


(19) 


(20) 


(21) 
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(4) Case 4, full cubic moment: 

C m (0) = C m< 6 + 3/2 C Ml 0 * + 2C Wj8 ^ !t (22) 

For cases 2, 3, and 4 it is required that C m (0) = 0. This 
actually does not constitute a restriction because the 
model reference axis can always be defined so that 0 — 0 
at the model trim position. 

The first phase of the analysis consists of the derivation 
of closed-form solutions for the dynamic stability deriva- 
tives of a vehicle exhibiting static aerodynamics as given 
by Eqs. (18-22). The first three pitching moment forms 
listed are completely analyzed. For the full-cubic moment, 
several case distinctions would be necessary for a com- 
plete analysis; only the physically most likely of these 
cases is considered. Nevertheless, the solution is much 
more complex than those corresponding to the other pitch- 
ing moments. This particular case, therefore, illustrates 
that for longer or more complicated functional forms of 
pitching moments an analytical approach is not practical, 
and numerical techniques should be employed. 

The limitation of closed-form solutions corresponding 
cnly to the pitching moments, given by Eqs. (19-22), 
imposes a serious restriction on the applicability of these 
solutions. However, the motion of a body exhibiting an 
arbitrary pitching moment can be analyzed by assuming 
one of the simpler moment shapes. The criteria and 
methodology for determining an approximating curve, 
given the original pitching moment and the free-flight tra- 
jectory, is the subject of the second phase of this analysis. 

As a verification of the validity and accuracy of the 
results presented here, the equations have been applied 
to two sample cases. A six-degree-of-freedom computer 
program was used to calculate trajectories corresponding 
to the two cases of nonlinear static aerodynamics. The 
equations derived in this report were used to calculate 
dynamic stability derivatives from the trajectory outputs 
of the program. Comparison of the calculated value of 
(C mfl + C m *) with the input value allows a direct verifi- 
cation of the accuracy of the formulas derived. 


II. The Energy Integral Relation and Its Application 

A. An Integral of Planar Motion 

The coefficients of the differential Eq. (17) are func- 
tions of 0 and 0\ and do not explicitly contain X. This 
suggests considering the motion in the phase plane (O' vs 
0 plane) rather than the coordinate plane (Fig, 4). This 




Fig. 4. Comparison of motion in coordinate 
and phase planes 


means eliminating dX in the differential equation in favor 
of terms involving dO Using the relations 0 " = dO'/dX 
and 0 ' = dO/dX and multiplying the differential equa- 
tion through by 9' dX, 


O'dO' 


pAd 


21 

X 


c m ( 0 ) do 


P Ad / pA\ dC m (0) 


pAd / pA \ 
21 \2m / 


dO 




+ | 5jC„(«) tt'di) 


pAd- 

2/ 


(C„, + C m .) e'de 

(23) 


For notational convenience, assume (without loss of gen- 
erality) the body has an initial angle — 9» < 0. Define 6 , 
as the positive angle to which the vehicle would oscillate 
in Vl> cycle if the system were completely conservative; 
i.e., if the static pitching moment were the only moment 
acting. Therefore, (9, is the angle that results in a net 
change of zero in potential energy resulting from C m (9) 
over Vs cycle, 

r C m (0)de~0 (24) 
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Equation (24) provides a method of calculating 0 U given the angle — 0 O and the pitching moment curve. Define BO 
as the loss in amplitude over V2 cycle as a result of the dissipated energy from the combined damping effects. In gen- 
eral, $!>> BO and, therefore, when integrating Eq. (25) over V2 cycle, the upper limit of integration ( 0 1 — BO) may be 
replaced by 0 X in all the second-order terms. 2 Also at the two amplitude peaks, — 0 O and (0 X — BO) the angular veloc- 
ity 0' = 0. Therefore, integrating Eq. (23) over V2 an oscillation cycle and solving algebraically for the dynamic sta- 
bility coefficient, 


(Cm q C w .) 


md- 

~ 


— md 

i J 

r“c.m»+'£r 

-0o J-e o 

o'* 

u 

de - [ Ol c D (0)0'd0 
J-Ot 1 

1 

f »> 

’ 0'd0 


-00 


(25) 


The equation for the damping coefficient in this form was first derived by Jaffe (Ref. 3) and has been called the energy- 
integral relation. 


To perform all the integrations indicated by Eq. (25) it is necessary to have an explicit expression for O' in terms of 0 . 
This can be obtained by ignoring all terms on the right side of Eq. (23) except the predominant pitching moment, 
as follows: 


O'dO' 


c m (o)de 


(26) 




(27) 


where k m ~ pAd/21. This estimate is an adequate representation in most physical cases. An iterative procedure to im- 
prove on it can be set up by using Eq. (27) as a first estimate, and repeatedly applying Eq. (23) to obtain progressively 
better approximations. However, the integrals become cumbersome very rapidly and, in general, the increase in accu- 
racy is insignificant and unnecessary. 

It will be convenient to state the energy integral equation in a standard form for reference purposes. Consider the 
first integral in the numerator of Eq. (25), as follows: 

tnd /m® 1-60) md F f 01 r&\- 6 e ~ | 

C m (0)d0 =-™[J Cm(0)d0 + J Cm(0)de J (28) 


The first term on the right vanishes by definition of 6,. The second integral can be evaluated approximately (to excel- 
lent accuracy) by using the trapezoid rule and the fact that 0, >> S0, 


md 

T/., 


c„, ( 0 ) de 


md , C„, (0,) + C m (0, — 86) 

— (-80) 2 


^ 0 C M 


(29) 


Equation (25) can be represented as follows: 




—j- S 0 C m ( 61 ) 


£ 


O'dO 


+ lift term — drag term 


(30) 


2 This is an approximation that does not need to be made when discussing numerical integrations later in the report. However, it makes no 
significant difference and, therefore, the point is not mentioned at that time. 
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For a wind-tunnel application, the first term on the right side usually comprises from 70 to 95% of the total. Therefore, 
the evaluation of the integral in the denominator of that term becomes the most important calculation when consider- 
ing the accuracy of the final answer. 

There is one additional equation that will be useful in succeeding sections. This is an estimate for the %-cycle dis- 
tance period of oscillation. By using 0 ' = d6/dX , 



Equation (27) may be used to substitute for 0' and the result integrated to obtain Xu 2 . This equation will be used in a 
later section to aid in the calculation of coefficients for an approximating moment curve from an experimentally mea- 
sured distance period. 

In a sense, Eqs. (25) and (27) comprise a complete solution to the general planar motion problem. Given the static 
aerodynamics in either functional or tabulated form and the angles ~0 O and (6h — 80) from an experimental oscil- 
latory history, these equations may always be integrated numerically to obtain a value for(C„, Q -f C OT .). However, it is 
still advantageous to solve the equations analytically for the most general case possible. A closed-form solution (in addi- 
tion to its ease of application) can be far more useful in parametric studies than a numerical integration and affords the 
investigator an immediate method to explore the relative importance of various terms in the differential equation of 
motion. The basic equations necessary for a complete analytical solution are Eqs. (24), (25), (27), and (31). The appli- 
cation of these equations collectively will be referred to as the energy-integral method. The next four sections of this 
report use the energy-integral method to obtain solutions for the dynamic derivative corresponding to four particular 
pitching moment curves. 


B. Offset Linear Pitching Moment: C MI (0) = C„, 0 + C mft 0 

I. Calculation of 0, from conservation of energy. 


By solving this quadratic for 0,, 


where 


1°' c m (6) do = o = c m< , 0 , + + c mo e„ 

0, = 0 n - 2 r 

r — C m JC mm 


(32) 


2. Estimate of 0'. 

O' = (2 *„)*[£ ( c «. + c m „e) do (2 *»)*[%* « 2 + C,„ a e ^)] ‘ 


(33) 


(34) 


3. Estimate of 'i- cycle distance period, X v 


X,.. 


(2 *:.)» 


do 1 r / (c„„s, + c.j \ 

+ c , + ( c „.„,. ' L"’ e c - - c -a > 2 J 


By substituting Eq. (33) for 0 , 


X,.. = 


(- k m C m Y a 


(35) 


(36) 
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4. Aerodynamic damping coefficient Using the lift and drag curves of Eq. (18) and the offset linear pitching 
moment and corresponding angular velocity in the energy-integral equation, the following solution results: 


. md J 289 md( C m 

(C mq + c ra .) — = - ^ ; _ r) — -j^- j + lift term - drag term 

Lift term = - w^r) + Ci « 0 + irh) - c ^ ( 2f - 3rL ' + 

+c, 2 [-y r* + J2 ~ r)2 + ^ 0n (° n ~ 2r) + ~ 0 6 -V M) ] 

Drag term = C Do - C Dl r + C„ 2 |\* + j(0 o - r) 2 J - C 0;( - -|e„r^l - ^„)J 

+ C„ 4 [f *- ( r + {(*« ~ r)') - H- r 2 0 o (r - { «. ) + *§ ( r - } fl. ) 2 ] 


(37) 


(38) 


(39) 


The first term on the right-hand side of Eq. (37) can be put into a simpler form by using Eq. (36) as follows: 

289 md / C>„ _ 4m 1 89 

"^«-f)T(T/ ~ (0„ — r) (40) 

5. Solution for axisymmetric shapes . For an axisymmetric body with its eg on the vehicle centerline, the drag curve 
would be symmetric about the line 9 — 0, and would, therefore, have an even-power series representation. The lift and 
pitching moment curves would be antisymmetric and have odd-power series expansions. Therefore, for an axisymmet- 
ric body with a linear pitching moment, Eqs= (37), (38), and (39) apply when C„, o , C L(} , C Ll , C Dl> and C D ^ are set equal 
to zero. The resulting formula is then identical to equations presented in Ref. 2. 


C. Bi-Linear Pitching Moment 

The most obvious, and perhaps the most useful, technique for dealing with a nonsymmetric moment curve is to 
approximate the plus and minus portions of the actual moment by two straight lines through the origin, as follows: 


C m (9) — C ma 9 , 0 — 0 

C m (0) = C ma+ 9, 0^0 


(41) 


Because in general, the slopes of the two curves will not be the same, this functional representation has a singularity at 
0=0. Therefore, all integral calculations must be performed separately in the two regions. 


I. Calculation of 9 X . 


/ o^ fei 

C m (9)d9= 0= / C ma _9 d9 + / 

«o J~ Q 0 «/° 


c m ode = -c„ 


■f + 


C~, -r- 


(42) 
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By defining t j- = 0j/05 = C m jC m the result is as follows: 


0 ! — ix6 u ~ 



( 43 ) 


2. Estimate of S'. Direct substitution in Eq. (27) yields the following two expressions: 

6’ = (-k m C m J*(6l - 0 2 ) 1 *, 0^0 

O' = (- k m Cm a J' A (81 - 8')* - ( - fc,„C m J* - 1 )*, 8 ^ 0 


(44) 


3. Estimate of the Vi-cycle-distance period, X % . 

l r de l 7 1 

X% = (-hnCnjuJ... (oi - 07* + ( - k m c m j»J„ (Or - 07* 


d0 


( — k m C Wa ) v - 2 ( — k m C nia J'A 2 ( k m C ma _Y A 


77/2 


j(1 + /*) 


(45) 


4. Aerodynamic damping coefficient. Evaluating the integrals of Eq. (25) in two pieces corresponding to 0 0 and 

0^0 and simplifying algebraically results in the following: 


mil- 


4m 1 80 


(C„, + C„,.) ^ -r=r- v + lift term - drag term 


2 C. 

Lift term f C L 

/At/ 0 


p/i #o 


(46) 


+ c,ot 


_f (i - * + /**) + b 


8 (l - 7(1 + 7) 




4 (i — u + «•-) 8 /ju - 1\ 

Drag term = C„ a 4- C« ( 0„ ^(/* - 1) + C c .,0? ^ + C„. ( 0 ‘ + ^ j 


+ 


C„ 4 0,’ 


(47) 


(48) 


D. Bi-Cubic Pitching Moment 

A natural extension to the idea of the bi-linear moment is to consider two approximating curves, each of which 
allows nonlinearities in its region of applicability. By using two distinct curves, rather than one higher order po y- 
nomial, a greater choice in the functional form is available. Choosing the two approximating curves as o -power 
series expansions results in a simplification of subsequent integrals. Therefore, this section deals with the fo owing 
moment curve: 


C,„ (0) = C m< 0 + 2 C m J\ 
C m (8) = C Ma 0 + 2C m J\ 


0—0 


(49) 
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1. Calculation of 0 X . 


f ‘ c m ( 0 ) dO = 0 - f° (C m 6 + 2C m 0 3 ) do + r (I Cm a O + 2 C m J*) dO 

J Oo J-*o Jo 

C„, J+ 0? + C m ^? - C.J - C m .Jl = 0 


or 


'■-{ 


- C„, tt ± [C ,4 + 4C„, 2+ (C„, a 05 + c„,._ 


2C„ 


F 


(50) 

(51) 

(52) 


The zh sign under the radical causes no ambiguity because the root desired is the smallest real positive root to Eq. (25). 

2. Estimate of O'. As in the previous case, the integral calculations must be broken into two parts corresponding to the 
two segments of the moment curve. For 0^0, 

9’ = (2*.)* [£* (CmJ + C m2 _0 3 ) dey = (kjc m j* (6* - OiY‘> («' + ^ + (53) 


For 0^0, 


0 ' = (2 k m )* | ~-c m J±- c,„., Y +f° (( c m e + 2 C m j 3 ) do y 


= (k„,yn- c u , a ( 9 5 - 0-) + c,„^0' - c„ a . 


By using Eq. (51) to simplify the latter equation, 

6’ = (k m C m J*(9 2 - 8l)»( 6* + ^ + AfV 

\ '-'Mi.,. / 


(54) 


(55) 


Equations (53) and (55) show that the angular velocity equation for the second l /4 cycle has the same functional form 
as the one for the first M cycle; they differ only by the replacement of 0 o by 0 Y and C m , z by C m2+ . This is as would be 
expected from the functional symmetry of the moment curve. Put another way, because the equation being used 
describes a conservative system, the kinetic energy at any given angle between — 0„ and 0 X must be the same, regard- 
less of whether the initial angle was — or 0 X . Therefore, the angular velocity at the given angle is also independent 
(except for sign) of which amplitude peak we start at. It follows from the symmetry in form of the two parts of the 
moment representation that the results from the first l A cycle will also apply to the second l A cycle of motion on sub- 
stitution for the pertinent constants. A similar argument holds for the other integrals involved in the calculations of 
the distance period and the damping coefficient. Therefore, for this moment, it is sufficient to consider only the first 
% cycle of motion, and to use these results with a substitution for the values of the integrals over the second cycle. 


3. Estimate of the l/ 2 -cycle distance period !, X % . Consider the first Vi-oscillation cycle, i.e., as 0 goes from — 0 o to 0, 


r — - - 

' J-e 0 O’ " “ 


= - (k m C m JK 


do 


(t 9 J - 0 


o) 1 '- ( 0 


+ ■ 




•r 


(56) 


This integral cannot be evaluated in terms of elementary functions. However, through use of a suitable substitution 
the integral may be transformed into the Legendre canonical form of an elliptic integral of the first kind (Ref. 8). A 
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case distinction and alternate substitutions must be made for the two possible algebraic signs of C m .,_, as follows: 

Case 1. Let > 0. The body is statically stable at — 0„, i.e., — C ma 8 0 — 2C„,„ O'f, >0. This implies the following three 
inequalities: 


2C 


+ 61 < 0 , 


-*■ + 61 < 0, 


w» /V 

- > 2 81 
HI., 


Let 

k - = 01/— (C ma /C m . t _ + 0%) '< 1, and <£ = sin 1 0/0 n to get the following: 


X, t = 



(57) 


where F (fc, </>) is the Legendre canonical form or normal elliptic integral of the first kind with modulus k and argument </>. 

Case 2. Let C m „ < 0, k~ - Ol/{C m JC hUl + 20S), </> = cos -1 0/0« 

The result is as follows: 

<58) 

Similar results can be obtained for the second \ \ cycle with C IH and 6 X substituted for C W>1 and 6 ih respectively. The 
‘•‘-cycle distance period is clearly the sum of the applicable l 4 -cycle results, X. 2 = X V| + Xt 4+ . Thus, there are four pos- 
sible expressions for the distance period, depending on the four possible sign combinations of the terms C m ., and C m;!+ . 
The results are summarized in Table 1. The necessary values of the elliptic integral F (k, tt/ 2) are readily available in 


Table 1 . Summary of distance periods over the Ye cycles 


First V 4 cycle 

Sign 

combinations 

ki 

x % - 

> o 
<0 

Q o 

> ( i r ( k 

c „ (t 

~C~ L + e ~° 

el 

(M' 4 V C„ , + C „ 2 / r \' , 2j 

'( \. V r(k -) 

C„, 

+ 2 $1 
m „ 

(M- \ c mtt +2c m3 el) V " 2 / 

Second V 4 cycle 

Sign 

combinations 

k; 

x. 4 * 

c „, (( > o 
c* it < o 

— 

e\ 

’ Y ’ V‘rb 

C m 

ci: • 

e\ 

c m 

-~ L + 2e\ 
1 

(k«)' 4 \ C,„ a +C m , #J/ \ 2/ 

(fc w )‘* V c m<t + 20 ^ 0 ;/ V ' 2 / 
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tabulated form. For reference purposes. Fig. 5 shows this 
function plotted vs the square of the elliptic modulus, k 2 . 


4 . Aerodynamic damping coefficient . Substitution of the 
bi-cubic moment and corresponding angular velocity in 
Eq. (25) again results in integrals that are nonelementary 
in nature. The most difficult integrations to perform are 
those which have the least influence on the accuracy of 
the solution, i.e., the integrals involved in the lift and drag 
terms. Furthermore, expressions for these two terms that 
are computed with the integrals assuming the bi-cubic 
moment are very cumbersome and difficult to apply. 
Therefore, because the loss in accuracy is small, the lift 
and drag terms will be handled in an approximate fashion. 
A succeeding section of this report derives methods of 
determining the coefficients of any particular (functional 
form) moment curve that will approximate an arbitrary 
moment curve. Therefore, given a bi-cubic moment, an 
approximating bi-linear moment may be determined and 
used to calculate the effects of lift and drag. Inspection of 
Eqs. (47) and (48) shows that the quantity 


P = = 9,/e 0 


describes the dependence on the pitching moment of the 
lift and drag terms. Because 6 0 and 6, are known, they 
may be used directly to calculate ft, and the approximating 
moment curve need not be actually displayed. (The physi- 
cal meaning and justification of this step will become 
clear in the following section.) Equations (47) and (48) 
will be used to approximate the lift and drag terms. 


The problem is, therefore, reduced to the evaluation of 
the denominator of Eq. (25), the integral. 



It is interesting, at this point, to investigate the magnitude 
of the error that would be introduced were a bi-linear 
moment used to approximate a bi-cubic moment in this 
remaining term of the energy integral equation. Figure 6 
(which also appears in Ref. 2) considers this question for 
a symmetric body, = C m , 2+ . The ordinate is the ratio 
of the prime term of the energy integral equation calcu- 
lated with equivalent cubic and linear moments, while the 
abscissa is the oscillation amplitude. The various curves 
represent different pitching moments from almost linear, 



Fig. 5. Complete elliptic integral of the first kind 
vs elliptic modulus 



0 20 40 60 80 100 


INITIAL AMPLITUDE 0 O , deg 

Fig. 6. Errors introduced by using a linear moment 
to approximate a cubic moment 

C m jCm 2 = —100, to fairly nonlinear, C m jC m . y = —5. As 
the figure shows, the error can become significant at the 
higher amplitudes. Although this plot concerns a sym- 
metric body, the magnitudes of the errors are indicative 
of the asymmetric case. 

Consider the integral, 



6'dO 


Again, there are four possible cases corresponding to the 
combinations of the algebraic signs of C„ K> and C nti>+ . An 
example with C m .,_ > 0 and C Wi , + < 0 will be carried out. 
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The other three cases may be generated by symmetry from the results given. 


1°' O'dO = (k m C m J* J" (0 2 - Ol)K + §^ + 0'sf 2 do + (k m C m J* £' (0 2 - 0\)* (V + §^ + Q'if do 


Let 

kt = - 

Then, 


01 


— + 01 


k 2 _ = 1 — k 'i , k ’i=- 


O'i 


fc; = 1 — fc;, </>. = sin 1 — 


+ 20? 


0o ’ 


( 59 ) 


= cos-' — 

Vy 


£ '* ■ {< c -- ) ' 4 [ - (fe + #i )] K [ (l + W ■ E (* ■ 1) - <• - « F (*- i)] 

+ (~C m J*(^ + 20?)* [(2*; - 1)E (ks£) + (1 - 


(60) 


where F (k, </>) and F ( k , <j>) are the Legendre canonical forms of elliptic integrals of the first and second kinds, respec- 
tively. This can be rewritten as in Eq. (61), letting Gj (kJ) and G, (k + ) have obvious denotations, as follows: 


£■ eve - g, (*.) + (c Mi+ )^ (^y Go (fc t )] 


(61) 


For reference purposes, Fig. 7 shows G, and G_ vs fc 2 . 

When the results are collected, the dynamic stability derivative is given by 

m( j- ~ (Cm a O\ + 2C„,_. t 0 1 ‘) 

(C,„ (( + C,„.) -y- = 7^7 + lift term - drag term 


[*' O'dO 


(62) 


where 



is given by Eq. (61), lift term by Eq. (47), and drag term by Eq. (48). 


E. Full-Cubic Pitching Moment 

The cases considered in the previous sections were intended to afford some degree of nonlinearity to the problem, 
and yet to limit the complexity of the integrals and resulting solutions. More general functional forms for the pitching 
moment, e.g., an unrestricted power series expansion, could be handled with the same approach. However, as the series 
becomes longer the solution becomes difficult and the number of needed case distinction multiplies. Therefore, direct 
numerical integrations of Eqs. (25) and (27) rather than an analytical approach, seem more practical. To exemplify the 
difficulties which can arise, consider the following pitching moment: 

C m (0) = C m< 0 + | C Wl 0 2 + 2C m .,Q'' (63) 
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Fig. 7. The functions Gi and G 2 
vs elliptic modulus 


This is the simplest power series representation avail- 
able beyond the pitching moments considered previously. 
A closed-form solution in terms of tabulated functions is 
obtained; however, the solution is restricted to one of sev- 
eral possible case distinctions and is more complex than is 
desirable. This particular pitching moment may be con- 
sidered as the longest power series expansion for which 
an analytical solution is practical. 


I. Calculation of 


J 0 ' c m (e)do= o = j (c„ a + c mi e, + c m j f) 


-T (c ma -c mi e 0 + c m2 et) 


or, rearranging 


V>m., \ / 

(64) 


One root of this equation must be — 0 O since 



C m (9) dO — 0 


Dividing by [6 X — ( — 0 o )j results in a cubic equation for 
0„ as follows: 


^- 9 > + fe-fe 9 * +9i > 

-(§^*.-§^ + «:)=° (« 5 ) 


An examination of the coefficients of this equation will 
reveal the following facts: (1) If C w , 2 is positive, there are 
three real roots, two positive and one negative; (2) if C,„ 2 
and C,„ ( are both negative, there is only one real root and 
it is positive; (3) if C m is negative and C, Wi is positive, 
there may be only one real root that is positive, or there 
may be three real roots, one positive and the other two 
either both positive or both negative. Over a large a range, 
the first case is the most likely physically; it is indicative 
of an initial increase in static stability with increasing «, 
and a subsequent leveling off and decrease in stability 
with further increase in angle. Because this is the most 
likely situation, the remainder of the discussion on this 
moment curve will be restricted to the first case. 


The algebraic solution of Eq. (65) using Cardan s formu- 
las (Ref. 9) is tedious and offers no simplification to the 
problem, as did the result for the bi-linear case. Therefore, 
expressions for the roots in terms of the coefficients will 
not be listed here. It is sufficient to say that the three roots 
of Eq. (65) may be obtained algebraically or by using a 
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simple iterative technique, such as the Newton-Raphson method (Ref. 10). For convenience, the Newton-Raphson 
method, as well as several other iterative techniques, are described in the Appendix. The oscillation amplitude peak 9 x 
will be the smaller of the two positive roots. The other two roots of Eq. (65) are useful in upcoming integrations; there- 
fore, assume this equation has been completely solved, either algebraically or numerically, for roots 0 U and — 0 A 
and that the following inequality holds: — 0 3 < — 0 O < 0 < 9 i < 0 2 . 


2 . Estimate of O'. By a direct application of Eq. (27), 

0' = (k m )Vt[c ma ( 0 * - 91 ) + C mi (0 3 4- e\ ;) + c m je* - 0;)]* 


( 66 ) 


In light of the prior discussion on roots of Eq. (65), this can be written in the more convenient form, 

0’ = (KCmJ* [(0 + 0.,) (9 + 0 O ) (0 - 0 X ) (9 - 0 2 )]tt 


(67) 


3. Estimate of ¥j-cycle distance period , X%. 

j_r 

c m ,)* [(0 + «.) (0 + 0„) (0 - 00 (0 - 0,)]' 4 


(k m C„ 


cld 


( 68 ) 


Let 


(0, ■+ 0o)(0, + 0a) . , f(0. + 0.0 (0 + 0(;)T- 

(0a + 0„) (0, + 00 ’ * ^ L (0, + 0„) (0 + 0a). 


Then 

X ' 4 [fcmC m2 (0 2 + 0„)(0 a + 0O]’'4 


(69) 


4. Aerodynamic damping coefficient. As in the bi-cubic case, the lift and drag contributions can be regarded as of sec- 
ondary importance and approximated with an equivalent bi-linear moment. Therefore, we are, again, left with the follow- 
ing term: 

P 0'do = P [(9 4- 0 8 ) (0 + 8o) {o - 9 1 ) (9 - 8 2 )]*d0 (70) 

J Of) 


It is at this point that the value of the closed-form solution becomes questionable. By using the substitution 


(0, -f 0Q(0 2 + 0Q 

(0 2 + 0») (0 1 + #0 


(0J + 0<o 

(0, + 0a) ’ 


0 < a- < k- 




r(0, + 0.0 (0 + 0.0 T' 

_( 0 , + 00(0 + 0 . 0 . 
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Elliptic integral parameter 
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the following solution results: 


/: 


0' do = ^°) (^* ^ ^°) (^ 3 ^°)~ /-. /% 7,2 


1(0* + Oo) (0 1 + 0,)]* 


G :i (a\k^) 


(71) 


where 


C, (a-, k 2 ) ~rj^ k 1 1 (17 > ® 2 > ^ ^ 1 (3fc 2 ~ — a 2 ) V ; , + (2a~k- + 2m- — 3k- — 77* ) V ; , + (a 2 — 1) (a 2 — k-) V 4 J 

V2 = 2W - iT(fc 2 - « 2 ) [“ 2£ ( fc> i) + ( * 2 “ 2) F ( fc> i) + (2 “* 2 + 2 “ 2 ~ 3k2 ~ “ 4) n (i ’**> fc )] 

V:< = 4(l-ffl)~(fr-S») [ k2F (*’ l) +2 G** + " 2 ~ 3fcL ’) 11 (| ^ *) + 3 & + 3fc2 “ - 2 « 2 ) V2 ] 


V 4 = 


6(1 -Tt- 


- _ — j^3fc 2 n , a 2 , k^ + 4(a-k- + a 2 - 3fc 2 ) V 2 + 5 (a< 4- 3fc 2 - 2 a 2 fc 2 - 2a 2 ) V, J 


and F (k, tt/2), E (fc, tt/ 2) and IT (7t/ 2, a 2 , fc) are complete elliptic integrals of the first, second, and third kinds, respectively. 
The dynamic stability derivative is, therefore, given by the following: 

nd, 80 (c ma 0, + | C mi Vi + 2 C„, 2 0? ) (0 X + 0 3 )% 

( Cm " + Cm ' a) ~T = 2 r KCZ (0, + W (07 + to) W< Z oj 2 G:< (a 2 , **) 


+ lift term — drag term 


(72) 


As can be seen, this is not an easy solution to apply to ex- 
perimental data. For convenience, the function G 3 (a 2 , k 2 ) 
is tabulated vs a z and k 2 in Table 2. Even with this aid, 
the solution (though workable) is not simple and repre- 
sents only one (the most likely one) of the possible case 
distinctions. Furthermore, a numerical approach has been 
suggested for obtaining the angles 0 l7 0 2 , and 0 3 . There- 
fore, it can be concluded that perhaps for this moment 
(and certainly for any pitching moment with a longer 
power series representation) either an approximating mo- 
ment curve or numerical integrations must be used. 


III. Equivalence Relations Between Moment Curves 

A. Equivalence Criteria 

The previous sections of this report have noted that 
there are a limited number of pitching moment represen- 
tations for which closed-form solutions are feasible. These 
include the offset linear, bi-linear, bi-cubic, and, possibly, 
the full-cubic moments that were discussed in Sections 
II-B to II-E. If we are given an arbitrary body, its pitching 
moment curve generally would not conform exactly to 


one of these representations. Therefore, if an analytical 
solution for the damping coefficient is desired, the actual 
pitching moment must be approximated with one of these 
four curves. This section deals with a method of deter- 
mining the coefficients of an approximating moment 
curve. 

Assume that an arbitrary pitching moment curve is 
given in either tabulated or functional form. An approxi- 
mating functional moment could be determined through 
use of conventional curve-fitting methods, such as the 
least-squares technique. However, with such an approach, 
although the approximating moment would be in some 
sense close to the original curve, there would be no assur- 
ance that any of the motion characteristics corresponding 
to these two curves actually would be the same. A better 
approach would be to use the matching of important pa- 
rameters of motion as the requirement for the selection of 
the coefficients of the approximating moment curve. The 
offset linear and bi-linear moments each have two coeffi- 
cients, and thus two characteristics of the angular motion 
can be matched. The bi-cubic and full-cubic moments 
each have three coefficients and allow the equating of 
three motion parameters. 
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The first (and most important) characteristic to match is 
the M>-cycle distance oscillation period, X Vi . In actual appli- 
cation, an average value of the decay per Vz cycle, 8(9, will 
be determined over several cycles of motion. It is, there- 
fore, necessary that the angular oscillations corresponding 
to the original and approximating moment curves remain 
in phase over this region of application. The second 
criterion is to match the amplitude peaks of the non- 
dissipative angular motion, i.e., given the initial angle 
~9 {h the corresponding positive peak 0^ should be the 
same for both moments. Matching of only these two con- 
ditions ensures great similarity between the motions, and 
will be the criteria used to determine the coefficients of an 
equivalent offset linear or bi-linear pitching moment. One 
more condition is available if either of the cubic moments 
is used on the approximating curve. Because the angular 
velocity 9' appears as a weighting factor of the dynamic 
stability derivative in the differential equation of motion, 
it would seem desirable to have similar angular velocity 
histories over the Vi cycle. This is borne out by the fact that 
in the energy integral equation, the most critical term is 

O' dO 

The second condition above ensures that the angular 
velocities will be identical (i.e., equals 0) at at least two 
points during the Vi cycle, — 0 n and at 6\. Furthermore, the 
matching of the Mi-cycle distance periods suggests that the 
angular velocity histories are of similar nature. The addi- 
tional degree of freedom afforded by the cubic moments 
will be used to match the maximum angular velocity over 
the M>-oscillation cycle, i.e., the angular velocity as the ve- 
hicle passes through zero angle-of-attack. These three 
conditions will then ensure a very good approximating 
representation of the original moment curve. 

The actual mechanics of determining the approximating 
moment curve depend on the knowledge of the true pitch- 
ing moment. It can be assumed that the original pitching 
moment is given either as a functional form, say a finite 
power series expansion, or as a tabulation. The polyno- 
mial representation may always be integrated directly, as 
indicated by Eq. (27), to yield O' as a function of 6 and 0„. 
By setting 9 — 0 in this expression, the maximum angular 
velocity, 0' max , as a function of 9 {) is obtained. An alterna- 
tive method to determine this quantity, which does not 
require a functional form for the original pitching moment, 
would be to numerically differentiate the experimental 
angular history at 9 = 0. However, the use of this tech- 
nique imposes more stringent sampling frequency and 


accuracy requirements on the experimental data. The re- 
sult of the integration mentioned can also be used to 
obtain 9 1 , the positive amplitude peak. At this point, the 
angular velocity will have to be zero. Then, 0 U is the small- 
est positive root to the equation 9 f — 0. This step can also 
be performed without actually knowing a functional rep- 
resentation for the pitching moment. A tabular form of 
the moment may be integrated mechanically or numer- 
ically, and the result used to determine 0 X . Finally, the 
^4-cycle distance period and the initial amplitude — 9 {) 
may be obtained directly from the free-flight trajectory. 
These four quantities, 0' max , ~0 o , 0 U and will be used 
to determine the approximating moment curve. Most of 
the necessary derivations have already been done in Sec- 
tion II. The following sections consist of manipulation of 
previous equations and statements of results. 

B. Offset Linear and Bi-Linear Pitching Moments 

1. Offset linear pitching moment: C m (9 ) = C„ l(( + C nin 0. 
Equations (33) and (36) provide the solutions for C„ I(1 and 

C, « rt , as follows: 


C„t a — 7 T~/k m (X^)~ 

(73) 

(VK 

(74) 

2. Bi-linear pitching moment: 


C m (0) = C ma J ; 0 — 0 


Cm (0) — Cm a+ 6; 0—0 


Rearrangement of Eqs. (43) and (45) yields the following: 


(75) 

c — c — 

(76) 


C. Bi-Cubic Pitching Moment 

C m (0) = CfnJ) + 2Cm 2 0 3 ; 9^0 

C m (0) = C ma O + 9^0 
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In the previous two cases, an actual calculation of the 
approximating moment curve was not necessary. The solu- 
tions for the damping coefficients, as given by Eqs. (37), 
(38), (39), (46), (47), and (48), are written in terms of the 
trajectory characteristics and do not explicitly involve 
the pitching moment. However, this is not the situation 
for the bi-cubic representation. The solution for the damp- 
ing derivative involves the pitching moment both explic- 
itly and implicitly in the moduli of the elliptic integrals. 
Therefore, it is necessary to produce the approximating 
moment curve. Furthermore, because of the implicit de- 
pendence through the elliptic integrals, the determining 
equations cannot be inverted directly to yield C m<t , C W2 , 
and C m2+ , but rather an iterative approach must be used. 
From Eq. (51) and from setting 0 ~ 0 in (53) or (54) to 
yield an expression for 0' max , 


equally valid, and it is the experimenters choice which to 
use. The inequalities which comprise this existence test are 
derived in the succeeding paragraphs. 

The simplest iterative procedure to use is different, de- 
pending on whether cubic coefficients are positive or 
whether at least one is negative. For a body with positive 
lift at zero angle-of-attack, it is likely that the restoring 
moment is less in magnitude at a given negative angle-of- 
attack than it is at the positive angle equal in absolute 
value. This means that C m . 2 _ is destabilizing, i.e., positive. 
The two cases with C W2 _ > 0 are discussed here. The re- 
maining two cases, C W2 < 0, could be handled similarly. 

I* C,» 2 > 0, Cm 2+ > 0. From Table 1 we get the fol- 
lowing: 


w - + c, "«\ e\ ) + e\ 


(77) 


f~> ifi max)* „ (0 max)* 

“ KOI “ Mr 


-C, 


(78) 


(U % [( c„ a + c m .ji ) F ( k ’ 2) 

+ ( _ Cm a + CZjf) F ( K l)] (80) 


where 


Instead of iterating on three equations in three unknowns, 
it is more expedient to solve the above two equations for 
two of the variables in terms of the remaining one, and 
then use the third equation to develop a recursion formula. 
Combining Eqs. (77) and (78) and rearranging 


k: - -C m . 2 Jl/{C ma + C m .Jl) 

k l = 

By using Eqs. (78) and (79), 


(C) ! (01 -e t) 01 
KOI 0* + 01 


(79) 



The third equation involves the ^-distance period X% and, 
therefore, will vary according to the algebraic signs of the 
cubic terms. A good estimate of the proper signs may be 
obtained by an inspection of the concavity characteristics 
of the moment curve. A precise method to determine the 
existence of a solution under the assumption of a particu- 
lar combination of signs is available through the two 
inequalities 0 < k 2 < 1 and 0 < k; < 1. The elliptic mod- 
uli, k 2 _ and k; may be written in terms of one unknown, say 
Cm 2 _, with the aid of Eqs. (78) and (79). Simplifying will 
result in inequalities imposing conditions relating the 
angles 0 o and 6 X . This is to be expected, i.e., the shape of the 
approximating moment curve would depend on the angles 
it must pass through when inscribing zero area. In general, 
more than one solution may exist; however, all should be 


where 

7 2 = Cm 2 - M 4 o 

(C ax) 2 

IO _ (01 -00 oioik m 

0? (0' ma *) 2 w - 

Examining the two inequalities implied by the elliptic 
moduli, the condition needed for this case to be applicable 
is as follows : 

4<-|<2 or 0.707 < j- < 1.414 
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The Newton-Raphson method is well suited for Eq. (81), for example: Let 


G(C m J X* 0' max F ( fc ’ 2 ) 0' inajI F (* + ’ 2) 


then 


dG 

dC,„.. 


*0 


dF 


H ) 


dk 


dF 


dk dC„ 


0 '„ 


2 ) 

dC m 


( 82 ) 


where 


dF ( k l) (1 - 

l(l-fc’) 




dk 


k m 0‘. 


d) K 


2 0-k 

0^1 


dC n 


2 (0' m «*) 2 k.’ dC,„. 2 _ 2(0' nmx ) 2 fc + 


The recursion formula is as follows: 


(Cm 2 _)i +1 (Gm 2 )i 


G(C ms J 


n 2 -' r dG(C m .,)/dc m ., 

Experience has shown that this iteration will, in general, converge rapidly, independent of the starting value. 
2. C„., > 0, C m ., + < 0. From Table 1 and the previous section, we get 

x,i = ( k ’ i) + (“ C7 a + 2C„, 2+ 0t) F ( fc+, 2)] 


(83) 


(84) 


where 




k m 0^ Crn 2 _ 

(«'m ax)> 5 


k: — 


C m 6i 


Cm„ + 2C m . 0 


By using Eqs. (78) and (79) 


Xl e r F (^ ’ 2) [■ 


met 


where k" is as above and 


k: 


(2 ei-ei)(6' max )>-k m o 


(oi-oiHQ'm Cm, 

(205 - ex){ 6 ' m „y-k m e\ei c m . 2 _ 


L'( fe ») 0 


(85) 
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Again, from examining the elliptic moduli, the condition for applicability is 0 < 0,/0 o < 1. Combining this with the pre- 
vious inequality, it follows that two solutions exist if 

l<^< 1.414 

Employing the Newton-Raphson method to solve Eq. (73) would require a derivative with respect to C wi2 , and would 
lead to a very complex recursion formula. A suggested approach is one which does not involve derivatives such as the 
method of accelerated iteration (see the Appendix) or perhaps a trial and error interval halving technique. Once a value 
for C„,., is obtained, it may be substituted in Eqs. (77) and (78) to calculate C Wft and C,„ 2+ . 

D. Cubic Pitching Moment: C m (0) = C m(x 0 + 3/2 C Wl 0 2 + 2C W2 0 3 

The solution for the cycle distance period again involves an elliptic integral and, therefore, an iterative approach 
is required. From Eqs. (64) and (66) we get the following: 


Cm. — 


1 U "«\0»+0»^ + 


r _ (fl'max) 2 . r . _ r fl2 

y»» a ~~ 7 /jo ' o Ml ( 7 0 o 


By using the latter two equations to solve for and C mi in terms of C nu> we obtain: 


C„i — 


~ Mi + ei) ~ Cro * M ' 


< 89) 

The distance period equation involves the auxiliary angles 0 2 and 0 3 . Two of the roots to Eq. (64) are known a priori, 
namely — 0 O and 0 X . A division then results in a quadratic for 0 2 and 0 3 > as follows: 

02 + r c mi -c ffl 2 (0o-0 1 ) -| $ + c ma - c mi ( 0 o - 0 Q + c W2 {oi - 0o0, + oi) = 0 (90) 

L Cm 2 J Cm 2 

By using Eqs. (88) and (89) to eliminate C nla and C w?1 from this quadratic, we get the following: 


b (£)’<*• - « 


$.1 — 
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The third equation is given by Eq. (69), as follows: 


Chi., 



1 

kn (02 + 0o) Wl 


+ 0 3 ) 


-0 


(93) 


responding to the initial angles — 0 € , were obtained from a 
plotted moment curve with the aid of a planimeter. This 
sample case illustrates the application of the analysis in 
the most general possible circumstances, i.e., when a func- 
tional representation of the pitching moment is not given. 


where 

(0i + 0 O ) (0, + 0 3 ) 

(02 + 0o) (0, +0a) 

Again, to avoid complicated derivatives, accelerated itera- 
tion or an interval halving approach can be used with 
Eqs. (92) and (93) to obtain C m 2 . C„i a and C m . x may then be 
calculated from Eqs. (88) and (89). 

IV. Applications and Accuracies 

A. Approach 

To exemplify the use of the formulas, and to verify their 
validity and accuracy, two sample cases have been worked 
out. The procedure employed was as follows: A six-degree- 
of-freedom computer program with aerodynamic coeffi- 
cients as tabular inputs was used to simulate a wind-tunnel 
free-flight trajectory. Each of the four pitching moment 
curves discussed in this report was used as an approxi- 
mating moment for analysis. Damping coefficients were 
calculated from the trajectory data corresponding to sev- 
eral initial amplitudes — 0 O . A comparison between the 
resulting values and the input value of (C niq + C W J is a 
direct measure of the accuracy of the equations. 

The aerodynamics for the two computed cases were 
chosen particularly to provide a severe test for the data 
reduction formulas. Both sets of static coefficients are very 
nonlinear and very asymmetrical. The curves are realistic 
because they are similar to proposed hypersonic lifting 
re-entry bodies; however, rather than being representative 
of such a class of configurations, they are more than typi- 
cally nonlinear in their static aerodynamics. 

For the first sample case, the static pitching moment is 
given by the cubic representation of Section II-E. The 
primary reason for using this exact representation was to 
provide a direct verification of the validity of the cubic 
analysis. The second sample case represents what may be 
considered as an arbitrary pitching moment curve, given 
only in tabulated form. The distance period and maximum 
angular velocity histories were obtained directly from the 
simulated trajectory. The positive envelope peaks 0, cor- 


B. Sample Case 1 

Consider a vehicle with static aerodynamic coefficients, 
as shown in Fig. 8, and an effective constant dynamic 
stability coefficient of — 1.0. The static pitching moment 
adheres to a cubic representation with C Wft — —0.1243, 
C Wl = -0.1471, and C m . t = 0.0678. Using the solution 
presented in Section II-E, the majority of the analysis may 
be done parametrically, vs initial angle —0 O> without 
actually using experimental data. This means that for a 
given pitching moment curve, all but the determination 
of, and multiplication by, 80 need be done only once, and 
can be completed prior to an experimental test program. 




ANGLE BETWEEN FREE -STREAM VELOCITY VECTOR 
AND MODEL CENTERLINE 0, deg 

Fig. 8. Static aerodynamics for sample case 1 
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< INITIAL AMPLITUDE 9 0 , deg 

Fig. 9. Roots of the potential energy equation 
for sample case 1 

Figure 9 shows the roots of Eq. (65), the potential energy 
equation, for this pitching moment. Shown is 0 1} the posi- 
tive amplitude peak corresponding to a negative peak 
— 0 O , 02 , and — 0 3 are the other auxiliary roots used in sub- 
sequent integrations. Figure 10 shows the two elliptic 
integral parameters, a 2 and k 2 , vs 0 O . These curves, in a 
sense, provide a measure of the nonlinearity of the pitch- 
ing moment. The larger the value of k 2 (i.e., the nearer to 
1) and the closer a 2 is to k 2 , the greater the nonlinearity of 
the problem. For this example, k 2 increases rapidly with 
0 O , and a 2 is always close to k 2 , therefore, indicating a very 
nonlinear pitching moment. 



Fig 10. Elliptic integral parameters for sample case 1 



INITIAL AMPLITUDE Q 0% deg 


Fig. 1 1 . Trajectory characteristics for sample case 1 

The maximum angular velocity 0' max and the %-cycle 
oscillation period X^ are plotted vs 0 O in Fig. 11. The 
dashed lines result from Eqs. (67) and (69) while the solid 
lines are from the computer-calculated trajectory. There is 
excellent agreement between the numerical and analytical 
solutions at all amplitudes. These trajectory characteristics 
are not actually needed to calculate the damping deriva- 
tive. They do, however, provide a standard of comparison 
for the various parts of the solution, and they become nec- 
essary later in the determination of approximate pitching 
moments. 



0 4 8 12 16 20 

INITIAL AMPLITUDE 0 O , d «Q 

Fig. 12. Integral of 0 ' for sample case 1 
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Figure 12 shows the integral of O' over V> of an oscilla- 
tion cycle. This term comes from Eq. (71) in conjunction 
with Table 2. Finally, Fig. 13 shows the lift and drag terms 
vs 0 O as calculated by Eqs. (47) and (48). Figures 8, 9, 12, 
and 13 may be used together with an experimental angular 
history to calculate dynamic stability derivatives. 




0 4 8 12 16 20 24 28 

INITIAL AMPLITUDE 6 0 , deg 


Fig. 1 3. Lift and drag terms for sample case 1 

Figure 14 shows a computer-calculated trajectory for 
this sample case. This figure points out a difficulty that may 
arise when dealing with a high-lift configuration. Because 
the lift is always positive, rather than of alternating sign as 
in the case of a symmetric body, the velocity in the Z 
direction increases steadily throughout the flight. Late in 
the flight (say during the sixth oscillation cycle), the 
angle (0 — «) has increased to about 1!1> deg because of 
this Z velocity. An angle of this magnitude after a Vi- 
oscillation cycle would not lead to difficulties because it 
has been adequately accounted for. However, considering 
the sixth cycle, itself, this term now enters as an initial con- 
dition (0 — a) 0 , which previously was assumed to be zero. 
The result is a smaller decay than would otherwise be 
expected when oscillating from a negative to a positive 
peak, and a larger decay when oscillating in the other di- 
rection. This, in itself, suggests a method to correct for 
this possible error. The decay in both oscillatory directions 
may be measured— an equivalent amount of decay at the 
positive peak corresponding to the observed decay at the 
negative peak determined from Fig. 9— and this equivalent 
decay averaged with the other observed decay to give an 



0 50 100 150 200 250 300 350 400 

MODEL POSITION RELATIVE TO MEDIA 
ALONG TUNNEL CENTERLINE X 


Fig. 14. Six-degree-of-freedom angular history, 
sample case 1 

effective value corresponding to an average 0„. This has 
been done for this particular trajectory, and the decay 
values were used in Eq. (72) to calculate damping coeffi- 
cients. The numerical results are presented in the first 
three columns of Table 3. The agreement with the input 
value of (C m + C m .) = —1.0 is excellent, within 1%, for 
all amplitudes considered. It must be noted, however, that 
the computer-determined decay, though rounded to the 
nearest 0.05 deg, is more accurate than could be deter- 
mined from a wind-tunnel test. In actual practice, a good 
procedure to follow would be to plot the positive and 
negative peaks, fair smooth curves through the plotted 
data, and determine the decay from the smooth curves. 
The decay could be plotted vs 0 O and smoothed again. 
With present data-gathering techniques, the result would 
probably be correct within 0.1 to 0.2 deg. This error in 
the determination of the decay would, of course, result 
in additional error in the final answer, the magnitude of 
which would depend on the actual level of the decay. 

For comparison purposes, this trajectory has been ana- 
lyzed assuming approximating pitching moments of the 
other three forms discussed in this report. For the offset 
linear and bi-linear moments, the damping derivatives can 
be calculated directly using Eqs. (37), (38), (39), (46), (47), 
and (48) without actually determining the approximating 
curve. The same is not true for the bi-cubic representation. 
Using the method outlined in Section III-C, the terms 
C M , C,„. , and C,„. (+ have been determined as functions of 
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INITIAL AMPLITUDE 0 O > deg 

Fig. 15. Equivalent bi-cubic pitching moment, 
sampie case 1 

<9 0 , and are shown in Fig. 15. Then, Eq. (62) was used to 
calculate (C m *• + The results for these three approxi- 
mating moments are listed in Columns 4, 5, and 6 of 
Table 3. The accuracy with the bi-cubic approximations 
remains fairly good, with the maximum deviation on the 
order of 10%. However, the two linear approximations fall 
off in accuracy rather badly, with the error approaching 
40% in the worst case. The sample case, then, empha- 
sizes the need for a nonlinear analysis when the pitching 
moment is indeed quite nonlinear. 




ANGLE BETWEEN FREE- STREAM VELOCITY VECTOR 
AND MODEL CENTERLINE 0, deg 

Fig. 16. Static aerodynamics for sample case 2 
C. Sample Case 2 

The static aerodynamics for Sample Case 2 are shown 
in Fig. 16. As in the first example, an effective constant 
dynamic stability derivative of — 1.0 was used. The static 
pitching moment is given only in tabulated or plotted 
form, i.e., without a functional representation. In the pre- 
vious example, the trajectory characteristics (9' ma x and 
could be generated from the known cubic moment and, 
therefore, a parametric analysis could be performed prior 
to experimentation. However, in this case, because an 
exact functional form for the moment is not available, it 
is necessary to have measured values of these variables 
before the determination of the approximating moment 


Table 3. Comparison of results using various pitching moment approximations 
for sample case 1 . Input value of (C m ^ + C m ^) is —1.0 


00, 

deg 

8$, 

deg 

(C„ + C„.) 

q a 

(cubic moment) 

«C. ? + C m .) 
(bi-cubic moment) 

<C„ + C„.) 

m q m a 

(bi-linear moment) 

(offset linear 
moment) 

10.0 

0.85 

— 1.011 

-1.076 

-0.926 

-0.860 

15.0 

1.25 

-1.004 

-1.091 

-0.871 

-0.751 

20.0 

1.50 

-1.002 

-1.110 

-0.778 

-0.624 
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curves. As determined from the computer-simulated tra- 
jectory, and 0' max > are plotted vs 0 Q in Fig. 17. The posi- 
tive envelope peaks corresponding to particular initial 
angles — 6 0 were obtained directly from Fig. 16 by using 
a planimeter, i.e., 

C m (0)d0 = 0 

The results are shown in Fig. 18. The amplitude decay 
history, determined from the computer calculations and 
Fig. 18, is shown in Fig. 19. Finally, the lift and drag terms, 




Fig. 17. Trajectory characteristics for sample case 2 


calculated according to Eqs. (46), (47), and (48), are 
plotted in Fig. 20. The group of Figs. 17-20 comprise the 
information necessary for the determination of the ap- 
proximating moment curves and the calculation of the 
dynamic stability derivative. 

The solutions corresponding to the two linear approxi- 
mating moments are written in terms of the trajectory 
variables and may be applied directly. The results are 
listed in Table 4, The offset linear moment does not pro- 
duce acceptable results; however, the bi-linear moment 



0 4 8 12 16 20 

INITIAL AMPLITUDE 0 O » d«Q 

Fig. 19. Amplitude decay history for sample case 2 



Fig. 18. Negative vs positive peaks for sample case 2 



ANGLE BETWEEN FREE- STREAM VELOCITY VECTOR 
AND MODEL CENTERLINE 0, deg 


Fig. 20. Lift and drag terms for sample case 2 
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Table 4. Comparison of results using various pitching moment approximations 
for sample case 2. Input value of [C m + C ms ) is —1.0 


00, 

deg 

80, 

deg 

< c „, + C m .l 

(cubic moment) 

(C % + C„.) 
(bi-cubic moment) 

< c„ g + C m .) 

(bi-linear moment) 

ic., + c..) 

(offset linear 
moment) 

5.0 

0.28 

— 1 .01 8 

-1.024 

-0.881 

-0.816 

10.0 

0.67 

— 1 .027 

-1.039 

-0.948 

-0.845 

15.0 

1.13 

- 1 .040 

-1.056 

-1,029 

-0.880 

20.0 

1.60 

— 1.046 

-1.075 

-1.056 

-0.867 



CO 



Fig. 21. Equivalent bi-cubic pitching moment, 
sample case 2 


appears to constitute a very adequate approximation for 
analysis purposes. 

By applying the iterative procedures described in Sec- 
tions III-C and III-D the two equivalent cubic moments 
were determined. These approximating curves are shown 
in Figs. 21 and 22. Equations (62) and (72) were used to 
calculate {C Mq + C OT<5 ). The results are listed in Table 4. 
The agreement with the input value of the damping coeffi- 
cient is quite good for both these moments, the maxi- 
mum deviation being less than 8%. 



CO 



Fig. 22. Equivalent cubic pitching moment, 
sample case 2 


V. Summary and Conclusions 

Over short time periods, dynamic stability derivatives 
generally have a relatively small influence on the free- 
flight angular motion of a vehicle. Observable effects con- 
sist of a slight convergence or divergence of the oscillation 
envelope. The deduction of a damping coefficient from a 
free-flight trajectory, therefore, requires a very accurate 
solution to the equation of angular motion. Employing 
assumptions characteristic to a wind-tunnel application, 
a solution for the dynamic damping coefficient, in terms 
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of integrals of the static aerodynamics and trajectory 
parameters corresponding to a nonaxisymmetric vehicle, 
was derived in this report. The format of the solution is 
such that the several quantities contributing to the vehicle 
amplitude change are completely uncoupled. Therefore, in 
the final evaluation of the integrals, approximations com- 
patible with the relative importance of each of these 
parameters may be employed in a term-wise fashion, 
thereby simplifying the solution without appreciable ac- 
curacy loss. In this manner, it was shown that the static 
pitching moment is the only aerodynamic coefficient in- 
volved in the primary term and, as such, constitutes the 
major limitation on nonlinearity of static aerodynamics. A 
systematic approach was applied to obtain closed-form 
solutions corresponding to several functional forms for 
the static pitching moment, each of increasing complexity 
and nonlinearity. The last pitching moment considered 
represents an upper bound with respect to functional gen- 
erality of the moment representation and practicality of 
obtaining and applying a solution. In addition, a method 
of determining an equivalent (for the purpose of analysis) 
moment curve conforming to any of the considered repre- 
sentations, given the free-flight trajectory and the vehicle s 
actual pitching moment (in either functional or tabulated 
form) was discussed. With this device, the equations de- 
rived herein may be applied to a body with completely 
arbitrary static aerodynamics. 

The examples demonstrate the validity of the results, 
as well as the computational procedure. The first example 
considers a pitching moment conforming exactly to a rep- 
resentation for which a closed-form solution is available. 
In such a case, a complete analysis can be performed para- 
metrically prior to actual experimentation. Furthermore, 


the most difficult portion of the computations, the deter- 
mination of equivalent moment curves, is not required. 
Calculations of dynamic stability derivatives from a simu- 
lated trajectory show excellent agreement with the input 
values, thereby verifying assumptions made in the devel- 
opment and solution of the equations. The second example 
considers a static pitching moment given only in tabulated 
form. In this instance, an a priori parametric analysis is 
not possible because trajectory information is necessary 
to determine an approximating moment curve. The choice 
of which functional moment form to use involves a trade- 
off. The two linear-type approximations require only very 
simple calculations, but do not lead to as good accuracy 
as do the nonlinear moments. On the other hand, the non- 
linear moments involve iterations and a digital computer 
for the computations. In some instances, this can lead 
to a great deal of time and effort in programming and 
applying an iterative scheme. All four approximating 
moments were determined and used to provide a com- 
parison of computational accuracies. In this particular 
example, the bi-linear approximation provides adequate 
results with a great deal less effort than required by the 
nonlinear moments and, therefore, is the logical selection 
for analysis. 

In a real application, it might be difficult to anticipate 
the accuracies of the various approximations solely from 
an inspection of the true static moment. A good procedure 
to follow would be to use a simulated trajectory to deter- 
mine the simplest approximating moment that could be 
used and still maintain the accuracy requirements. In any 
instance, it is expected that at least one of the moment 
forms and corresponding solutions will yield results accu- 
rate to within 10%. 
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Appendix 

Iterative Formulas for the Solution of f (X) = 0 


1. Method of false position. 


. _ xrm 
i+l m-fw 


4. Method of accelerated iteration. To apply this 
method, it must be possible to rewrite the equation in the 
form f (X) = g (X) — X = 0. 


2. Newton-Raphson method. 


Let 


X i+1 =X* 


/(Xi) 

f(x i) 


Xij = g (Xi), X 2j = g (X, .) 


3. Newton's abbreviated method. 


Then, 


• /(Xi) 

Y(o) 


_ x<x 2j — (X ai ) 2 
Xi + X 2; - 2X 1; 


Nomenclature 


A 

model reference area, 7rd 2 /4 

Cy 

normal force coefficient, normal force/ 

a 2 

elliptic integral parameter 


q*A 

c A 

axial force coefficient, axial force/q^A 

d 

model diameter, reference length 

C D 

drag coefficient, drag fore e/q^A; C Dq , 

D 

drag force 


C Di , C, v C D:<y c» 4 are coefficients in 
power series expansion for drag co- 
efficient 

&7> 

E(k,<t>) 

center of pressure offset from reference 
axis 


elliptic integral of the second type 

c. 

lift coefficient, lift fore e/q^A; C L{)> C La , 
C Ll> C/ 2 are coefficients in power series 

F (k, <t>) 

elliptic integral of the first type 


expansion for lift coefficient. 

F, 

axial force 


static pitching moment coefficient, pitch- 
ing moment/ q^Ad; C ni{} , C ma , , C mv Cm a 

F.v 

normal force 

acceleration due to gravity 


are coefficients in power series expansion 

g 


for static pitching moment coefficient. 

I 

model moment of inertia about a trans- 

(C«, + C m .) 

dynamic stability coefficient; 


verse axis at center of gravity 


dCm/d (^Pj + dCm/d 

k._ 

modulus of elliptic integrals; k ± is equal 
to complementary modulus, k 1 — 1 — k 2 
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Nomenclature (contd) 


pAd/2I 

l ( .y distance from model nose to center of 
gravity 

l r{) distance from model nose to center of 
pressure 

L lift force 

m model mass 

M rs sum of moments about the center of 
gravity 

q w free stream dynamic pressure 
t time 

V model velocity relative to media 

Voo free stream velocity 

x m model position relative to inertial system 
along tunnel centerline 

X model position relative to media along 
tunnel centerline 


Z 9 z 


m 


so 

e 

On 

0 , 




? 

</> 

(•) 

(') 


model position in vertical direction 
model angle-of-attack 
amplitude decay 

angle between free-stream velocity vec- 
tor and model centerline 

initial amplitude 

amplitude of conservative system after 
V 2 cycle 

elliptic integral of the third type 
gas density 

amplitude of elliptic integrals 
derivative with respect to time 
derivative with respect to distance 
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